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ABSTRACT 


This paper describes an improved implementation of the 
Chebyshev ephemeris generation capability in the opera- 
tional version of the Goddard Trajectory Determination 
System (GTDS) . Preliminary results of an evaluation of 
this orbit propagation method for three satellites of 
widely different orbit eccentricities are also discussed 
in terms of accuracy and computing efficiency with respect 
to the Cowell integration method. An empirical formula is 
also deduced for determining an optimal fitting span which 
would give reasonable accuracy in the ephemeris with a 
reasonable consumption of computing resources. 
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SECTION 1 - INTRODUCTION 


This document presents an improved implementation of the 
Chebyshev ephemeris generation capability in the Goddard 
Trajectory Determination System (GTDS) . The reimplementa- 
tion was necessary to resolve a System Failure Report on 
the operational version of GTDS and to improve the clarity 
of the computer program code to make it more readable and 
maintainable. The improved implementation employs the 
same Chebyshev polynomial/Picard iteration scheme as pre- 
viously implemented (described in References 1 and 2) but 
exhibits a marked improvement in accuracy and efficiency 
(see Appendix B) . The improved implementation fits the 
Chebyshev polynomial to satellite ephemeris data displaced 
as a function of time in accordance with the roots of the 
Chebyshev polynomial. This displacement is dependent on 
the degree of the polynomial. 

The advantages of using Chebyshev polynomials as inter- 
polating polynomials and the computational scheme in GTDS 
are briefly described in Section 2. Section 3 discusses 
general application of the improved implementation of the 
Chebyshev method to orbits over a wide range of eccentric- 
ity. The results are analyzed in deducing an empirical 
formula for determining an optimal fitting span that would 
consume a reasonable amount of computer resources and 
still provide a reasonably accurate ephemeris. A brief 
summary of conclusions is presented in Section 4. 

Appendix A briefly discusses the properties of Chebyshev 
polynomials, the formulation of an interpolating poly- 
nomial consisting of a linear combination of Chebyshev 
polynomials of different degrees to represent accelera- 
tion, and the integration of the interpolating 
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polynomial to generate satellite ephemerides. Appendix B 
contains the results of a comparison of the new and pre- 
vious implementations of the Chebyshev ephemeris genera- 
tion method in GTDS. 
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SECTION 2 - THE USE OF CHEBYSHEV POLYNOMIALS 
TO GENERATE EPHEMERIDES 


2. 1 ADVANTAGES OF USING CHEBYSHEV POLYNOMIALS AS 
INTERPOLATING POLYNOMIALS 

In principle, any function characterized by a discrete set 
of values can be approximated by a polynomial or a linear 
combination of polynomials. Such polynomials may be ex- 
pressed as Chebyshev polynomials, Legendre polynomials, 
Laguerre polynomials, or any other polynomial form which 
is expedient for mathematical/computational analysis. For 
instance, in the case of a satellite trajectory, the posi- 
tions at a series of selected times determine a polynomial 
consisting of a Chebyshev series within the time inter- 
val. The significant advantages of using Chebyshev poly- 
nomials to fit a satellite trajectory are that the error 
in the approximation is distributed evenly over the inter- 
val and that the maximum error is reduced to the minimum 
or near-minimum value (References 3 and 4) . 

Once this interpolating polynomial is established, the 
position of any other time within the interval can oe 
easily interpolated. If a long ephemeris is to be stored 
for any reason, it is plausible to use a small amount of 
computer storage to store only coefficients for the inter- 
polating polynomial instead of using a large amount of 
space to store the entire ephemeris. One familiar example 
is the Solar/Lunar/Planetary Ephemeris File (SLP File) , 
which is stored as coefficients of Chebyshev polynomials 
for GTDS and other trajectory determination systems to 
interpolate noncentral body positions for evaluating per- 
turbations on a satellite. Another possible application 
would be to store the coefficients of Chebyshev poly- 
nomials to represent the ephemeris of a Tracking and 
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Data Relay Satellite (TDRS ) in the onboard computer of a 
user satellite for autonomous orbit determination. 

2.2 COMPUTATION SCHEME IN GTDS 

In order to apply the mathematical theory described in 
Appendix A, one must know the acceleration, x(£), as a 
function of time to fit a Chebyshev interpolating poly- 
nomial. However, this is not the case for near-Earth 
spacecraft because of the nonlinearity of the perturbing 
forces, namely x depends on x which is in turn determined 
from x. Therefore, the Picard iteration method is used in 
GTDS to incorporate the Chebyshev series ephemeris genera- 
tion method. The computational procedure is described in 
the following paragraphs. For discussions related to the 
mathematical aspects of the method, see Reference 1. 

Suppose an ephemeris is requested from t to t with a 

d Z 

fitting span (or equivalent step size) of H which is equal 
to (t^ - t ). The entire ephemeris will consist of a 
series of spans which are represented by different 
Chebysnev interpolating polynomials. The default fitting 
span in GTDS is 5400 seconds. The allowable range of the 
degree of the Chebyshev interpolating polynomial is from 4 
to 48 with a default of 36. 

Within a fitting span, the roots (£ k of the Chebyshev 
polynomial of the highest degree plus 1, (n + 1) ) in the 
interpolating polynomial are first computed according to 
Equation (A-7) . These roots are then transformed back 
into time, i.e., £ k + t k , k = 1, 2, ..., n + 1. 

GTDS uses boundary conditions at the beginning of the fit- 
ting span, i.e., the position and velocity at t to 
obtain positions and velocities at t-^, t 2 , . .., t k , 

. t n+ i with a two-body central force field to start 
the iteration scheme. With the positions and velocities 
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at available, the perturbations, x(^ k ) can now be 
estimated at these instants and the Chebyshev coeffi- 
cients, C., are subsequently computed using Equa- 
tion (A-ll) . At this point, the Chebyshev interpolating 
polynomial for the acceleration. Equation (A-8) , is es- 
tablished. 

The next step is to successively integrate the Chebyshev 
interpolating polynomial twice according to Equa- 
tions (A-14) and (A-17) to obtain interpolating poly- 
nomials, Q n+1 and R n+ 2 * for velocity and position, 
respectively, in the fitting span (t a , t^) . The posi- 
tion and velocity with perturbations included at the end 
of the fitting span, t^, or any other time can be easily 
interpolated. The first loop of the iterative scheme is 
essentially completed at this point. 

In the next loop, GTDS uses positions and velocities in- 
terpolated from the interpolating polynomials, Q ^ and 
R n+2 f at the roots to estimate acceleration. After fit- 
ting the polynomial to the accelerations, it is again in- 
tegrated twice to obtain polynomials for velocity and 
position. The position interpolated at the end of the 
fitting span in this loop is compared with that obtained 
in the previous loop. 

This iterative scheme is repeated until the differences of 
the position components of the two successive loops at 
t^ are less than a tolerance (default value = 

10 -6 kilometers) . At this moment, the fitting procedure 
for the span ( t , t b ) is completed. 

After ephemerides are generated and the Chebyshev coeffi- 
cients for velocity and position are optionally saved, the 
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fitting span is advanced one step forward to 

(t b , t b + H) . This scheme is continued until all the 

spans are fitted. 


SECTION 3 ~ APPLICATIONS OF THE IMPROVED IMPLEMENTATION OF 
THE CHEBYSHEV EPHEMERIS GENERATION METHOD 


The improved implementation of the Chebyshev ephemeris 
generation method is applied to satellites of different 
orbital eccentricities to study the behavior of the 
Chebyshev polynomial representation in order to find an 
optimal set of parameters, such as fitting span and degree 
of the Chebyshev polynomial, for different satellites. A 
series of computer runs on GTDS with the new Chebyshev 
implementation was obtained. The ephemerides from the 
Chebyshev ephemeris generation method are compared with 
those from the Cowell integration method in terms of ac- 
curacy and efficiency. The results are discussed 
separately for a near-circular orbit, an elliptical orbit, 
and a highly eccentric orbit in the following sections. 

An attempt to find an empirical formula for determining 
the optimal fitting span for these orbits is also dis- 
cussed. 

-3 

3. 1 NEAR CIRCULAR ORBIT (ECCENTRICITY = 10 ) 

The GEOS-3 satellite was chosen for this case study. The 
eccentricity of the GEOS-3 orbit is 0.00098 and the semi- 
major axis is 7225 kilometers. The fitting spans used in 
this case range from P/4 to 2P, where P is the period of 
the satellite. For each fitting span, several runs with 
different degrees of Chebyshev polynomials were made. The 
ephemeris of every run was compared by using the GTDS 
Ephemeris Comparison Program with the reference ephemeris 
generated by the Cowell integration method with a 24-sec- 
ond step size using perturbations identical to those used 
in the Chebyshev method. The maximum differences in posi- 
tion vector, I A R | , between the two 
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ephemerides are plotted in Figure 3-1 as a function of the 
degree of Chebyshev polynomials and the fitting span. 

The maximum difference decreases very rapidly as the de- 
gree of Chebyshev polynomials increases. However, after 
reaching a critical degree of the Chebyshev polynomials, 
the maximum difference bottoms out and does not decrease 
any further. 

The saturation of the maximum difference occurs at a lower 
degree of the Chebyshev polynomials for a shorter fitting 
span. This saturation level generally increases with the 
fitting span. 

Since the step size of -24 seconds used in the Cowell inte- 
gration method in generating the reference ephemeris is 
relatively very small, the maximum difference in position 
vectors between the Chebyshev and Cowell ephemerides can 
be loosely regarded as the accuracy of the fit of the 
Chebyshev polynomials. Therefore, Figure 3-1 demonstrates 
one significant phenomenon: once the saturation level is 

reached, for a particular fitting span, adding higher de- 
grees of the Chebyshev polynomials not only does not im- 
prove its accuracy, but decreases its efficiency. This is 
further evaluated by examining the computer resources, 
mainly CPU time, consumed by each of the computer runs. 

All the runs for GEOS-3 satellite were executed on the 
GSFC IBM S/360-75 Cl computer. However, some of the runs 
were executed in the "low-speed" core of the CPU, which is 
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roughly three times slower than "high-speed" core. Also, 
there are a number of runs executed partly in low-speed 
core and partly in high-speed core. The CPU time consumed 
depends on the proportion of high-speed core and low-speed 
core used. This nonuniform CPU scale- makes the comparison 
not so straightforward. Instead of reexecuting these runs 
in high-speed core, the CPU time of these runs are cali- 
brated through force model calls, as described below. 

In GTDS , the "Number of Times Forces Called For Full 
Model" in the statistics report is provided at the end of 
a run. For a numerical integration method, such as the 
Cowell method or the Chebyshev method, the full perturbing 
force, including harmonic geopotential field, noncentral 
body gravitational field, and nonconservative forces, is 
evaluated at each integration grid point according to the 
options specified. The number of times the full perturb- 
ing force is evaluated is proportional to the CPU time 
used in a run. In Figure 3-2 this number is plotted 
against CPU time for only those runs executed in high- 
speed core. Although the points plotted are somewhat 
scattered, there is a linear relationship between this 
number and the CPU time. 

For comparison, the number of times the full force is 
evaluated (7236 times) and the CPU time (0.85 minute) are 
also plotted in Figure 3-2 for the reference run of Cowell 
method with a 24-second step size. It is interesting to 
note that the Chebyshev method with any reasonable accu- 
racy is much slower than the Cowell method. Consequently, 
for ordinary purposes other than those that require 
Chebyshev coefficients, it is at least not recommended to 
use the Chebyshev method to generate ephemerides for a 
spacecraft of circular orbit at a lower altitude. 
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NUMBER OF TIMES FULL FORCE IS EVALUATED 


"EXECUTED ON THE IBM S/360-75 Cl COMPUTER 


Figure 3-2. Calibration Curve of CPU Time Consumed 
by Computer Runs Used for This Study 
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After the CPU time of the runs which were executed partly 
in high-speed core are calibrated by using this linear 
relationship, ^ the CPU time of all the runs of Chebyshev 
method is plotted against the degree of the Chebyshev 
polynomials in Figure 3-3 for different fitting spans. 

The CPU time consumed is approximately linearly propor- 
tional to the degree of the Chebyshev polynomials. The 
CPU time is also plotted in Figure 3-4 against the fitting 
spans for degrees 18, 28, 38, and 48. 

The curves in both Figures 3-3 and 3-4 give the impression 
that the fitting span of one satellite period would be the 
most desirable one to use for the Chebyshev method as far 
as CPU time is concerned. However, the accuracy of the 
fit may not be desirable for the situation. For this rea- 
son, the accuracy information is also included in Fig- 
ures 3-3 and 3-4 by different shadings of the plot symbols 
to avoid the possibility of drawing misleading conclusions. 

Since the accuracy of the Chebyshev method bottoms out at 
a critical degree of the polynomial (Figure 3-1) and the 
CPU time used increases linearly with the degrees of the 
polynomials, a trade-off can be performed to study the 
benefit or penalty of using a higher degree than is neces- 
sary . 

The results of the trade-off study are presented in Fig- 
ures 3-5, 3-6, 3-7, and 3-8 for fitting spans P/4, P/2, P, 
and 2P, respectively, by combining the results in Fig- 
ures 3-1 and 3-3. The CPU time or the accuracy is normal- 
ized with respect to that of a data point 


1 This calibration curve is not necessarily valid for other 
experimental conditions or other satellites. 
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Figure 3-3. CPU Time Versus Degree of Chebyshev Polynomials 
for Various Fitting Spans for GEOS-3 
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Figure 3-7. Changes in Relative Accuracy and Efficiency for 
a Fitting Span of P for the GEOS-3 Orbit 
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Changes in Relative Accuracy and Efficiency for 
a Fitting Span of 2P for the GEOS-3 Orbit 
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corresponding to the critical degree of the polynomial on 
the saturated portion of the accuracy curve. For example, 
the CPU time used in a computer run for fitting a 
Chebyshev polynomial of the ith degree over a span of P/4 
is normalized with respect to the CPU time consumed for 
fitting 13th-degree Chebyshev polynomials over the same 
span, i.e.. 


[(cpu). - (cpu) 13 1 / [ (cpu) 13 ] 

Likewise, the accuracy of the fit is also normalized with 
respect to the 13th-degree Chebyshev polynomials, i.e., 

^| A ^lmax,i " I A ^lmax,13^ | A ^|max,13^ 

For the fitting span of P/4, CPU usage doubles without any 
benefit at all when the degree is increased from 13th to 
25th. Actually, the accuracy has deteriorated by about 
10 percent. If the degree is reduced from 13th to 12th, 
the CPU consumption saved is only 0.6 percent, but the 
penalty is a significant 70 percent decrease in accuracy. 
Therefore, it is very desirable to predefine the require- 
ment for support to be accuracy-bound or CPU-bound for 
selecting the fitting span and the degree of the Chebyshev 
polynomials. An arbitrary combination of these parameters 
may either produce an ephemeris with accuracy so poor that 
it is not usable or consume more computer resources than 
necessary. 

Another area of trade-off consideration is whether support 
is accuracy-bound or storage-bound. The total number of 
Chebyshev coefficients is directly proportional to the 
degree of the Chebyshev polynomials and inversely 
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proportional to the fitting span over a predefined arc 
length. If these coefficients are to be saved in a 
limited amount of space for general applications, such as 
ephemeris representation on an onboard computer for satel- 
lite navigation or autonomous spacecraft, an appropriate 
combination of degree and fitting span must be selected 
for an efficient usage of the storage within a required 
accuracy constraint. 

3.2 ELLIPTICAL ORBIT (ECCENTRICITY = 0.1) 

The IMP-7 spacecraft, with an orbit eccentricity of 0.11 
and a semimajor axis of 223,670 kilometers, was selected 
to represent the elliptical orbit. Three sets of computer 
runs were obtained for fitting spans of P/4.5, P/2, and 
P. The results are shown in Figure 3-9. 

The behavior in the variation of accuracy with the degree 
of the Chebyshev polynomials is essentially the same as 
that shown in Figure 3-1 for a near-circular orbit. The 
accuracy improves very rapidly as the degree increases and 
then saturates after a critical degree is reached. 

3.3 HIGHLY ECCENTRIC ORBIT (ECCENTRICITY = 0.9) 

The ISEE-1 spacecraft orbit, with an eccentricity of 0.91 
and a semimajor axis of 75,500 kilometers, was selected as 
representative of a highly eccentric orbit. With a fit- 
ting span of P/4, equivalent to 51,600 seconds, the best 
accuracy of the ephemeris represented by Chebyshev poly- 
nomials of the 48th degree over one satellite revolution 
(equivalent to a 2.4-day arc) is 228 kilometers with re- 
spect to the reference ephemeris generated by the Cowell 
method. An ephemeris with accuracy this poor may not be 
very useful. 
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3-9. IMP-7 Ephemeris Accuracy for the Chebyshev Poly 
nomial Fit Compared With the Cowell Method 



Further tests were conducted with drastically reduced fit- 
ting spans of P/40 and P/80. The results are presented in 
Figure 3-10. The accuracy obtained was comparable to that 
shown in Sections 3.1 and 3.2. As in the cases of cir- 
cular orbit and elliptical orbit, the- accuracy curves show 
the similar behavior in the variation of accuracy with the 
degree of the Chebyshev polynomials. 

3.4 AN EMPIRICAL FORMULA TO DETERMINE THE OPTIMAL FITTING 
SPAN 

Ideally, in applying the Chebyshev method, one would like 
to obtain the highest accuracy with a minimum amount of 
CPU time for the lowest possible degree and the longest 
possible fitting span. However, so straightforward an 
application is not possible because those factors compete 
with each other in a rather complicated fashion as demon- 
strated in Sections 3.1, 3.2, and 3.3. An attempt was 
made to find an empirical formula for determining an opti- 
mal fitting span in terms of satellite period. 

From Figures 3-1, 3-9, and 3-10, it is obvious that a 
longer fitting span requires a higher degree for the 
Chebyshev polynomials in order to achieve acceptable fit- 
ting accuracy, i.e., the fitting span should be propor- 
tional to the degree of the Chebyshev polynomials. Fur- 
thermore, the fitting span must be substantially smaller 
for a highly eccentric orbit than for a circular orbit. 
From these arguments, a very crude empirical formula re- 
sults: 


H 


C 


DP (1 - e) 


2 


(3-1) 


where H = the fitting span of Chebyshev polynomials in 
terms of satellite period 
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Figure 3-10. ISEE-1 Ephemeris Accuracy for the Chebyshev 

Polynomial Fit Compared With the Cowell Method 
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D = the degree of the Chebyshev polynomials 
P = the satellite period 
e = the eccentricity of the orbit 
C = an empirical constant 

The value of the constant, C, depends on the degree of the 
Chebyshev polynomials, 


C 


1 

40 


D s 20 


C = 


1 

20 


D > 20 


To demonstrate the validity of this empirical formula, the 
following examples are given and results are shown in 
Table 3-1. 

To represent the GEOS-3 ephemeris (e 0) with Chebyshev 
polynomials of the 10th degree, the fitting span computed 
using Equation (3-1) is 2P, which gives an accuracy of 
0.27 meter over 28 periods (2 days). If Chebyshev poly- 
nomials of the 20th degree are desired for GEOS-3, the 
fitting span given by the empirical formula is P/2, which 
gives an accuracy of 0.13 meter. 

For IMP-7 (e 0.1), the fitting span computed from the 
empirical formula for a 40th-degree Chebyshev polynomial 
is roughly 3P/2 with an accuracy of 20 meters. For a 
20th-degree Chebyshev polynomial, the fitting span would 
be P/2.5, giving an accuracy better than 0.2 meter. 

In the case of ISEE-1 (e 2 0.9), fitting a 48th-degree 
Chebyshev polynomial requires a fitting span of P/40 to 
achieve 0.02 meter accuracy. For a 30th-degree Chebyshev 
polynomial, a fitting span of P/70 gives an accuracy of 
0.03 meter. 
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Table 3-1. Ephemeris Accuracies for Fitting Spans Determined Using 
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The empirical formula is applied to the Tracking and Data 
Relay Satellite (TDRS) and the results are also included 
in Table 3-1. The TDRS is to be a geosynchronous satel- 
lite with an eccentricity of nearly zero and a semimajor 
axis of 42,000 kilometers. The fitting span computed from 
the empirical formula for a 20th-degree Chebyshev poly- 
nomial is P/2 which gives an accuracy of 0.14 meter for 
the ephemeris over a 31-day arc of 31 revolutions. If a 
40th-degree Chebyshev polynomial is chosen, the computed 
fitting span is 2P which gives an accuracy of 0.16 meter 
over the same arc length. 

To further verify the validity of the empirical formula, a 
nonexistent Satellite-X with an eccentricity of 0.5 and a 
semimajor axis of 13,200 kilometers was tested. The peri- 
gee height is 6,600 kilometers, about 200 kilometers above 
the surface of the Earth, and the apogee height is 
19,800 kilometers. The relative importance of all per- 
turbing forces, such as a higher-order harmonic geopoten- 
tial field, atmospheric drag, and solar radiation 
pressure, exerted on Satellite-X varies at different 
positions on the orbit causing the magnitude of the 
trajectory variation to differ along the orbit. Near the 
perigee, a shorter fitting span and a higher degree of 
Chebyshev polynomials may be needed to meet required 
accuracy criteria because of the effects of a higher-order 
harmonic geopotential field and the atmospheric drag. 

Near the apogee, a medium fitting span and medium degrees 
of the Chebyshev polynomials may be required because of 
the large curvature of the trajectory in combination with 
the trajectory variation due mainly to the solar radiation 
pressure. While in the vicinities of 90 degrees and 
270 degrees of anomaly of the orbit, the trajectory is 
rather linear and a lower degree and a longer fitting span 
may be sufficient. It is not possible, however, to apply 
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several different fitting spans and degrees over each 
revolution in a single computer run setup with the current 
GTDS, which only allows a uniform fitting span and a sin- 
gle choice of degree for the Chebychev polynomials. 

Two test runs were made with fitting spans of P/8 and P/2 
computed by using the empirical formula. Equation (3-1) , 
for Chebyshev polynomials of the 20th and 40th degrees, 
respectively. The accuracy for a P/8 fitting span with a 
20th-degree Chebyshev polynomial is 0.55 meter, and 
0.62 meter for a P/2 fitting span with a 40th-degree poly- 
nomial over a two-day arc of 11.5 revolutions. 

With the exception of the case of the 40th-degree poly- 
nomial with a 3P/2 fitting for IMP-7, all the cases seem 
to favorably support the validity of the empirical for- 
mula. However, the formula still should be used with ex- 
treme caution, perhaps only as, a rough guideline to 
establish a preliminary set of parameters for the 
Chebyshev ephemeris generation method. 
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SECTION 4 - CONCLUSIONS 


The Chebyshev ephemeris generation method is reimplemented 
in the operational version of GTDS . The conclusions from 
the testing results for this new implementation are sum- 
marized below. 

• The new implementation is more efficient and pro- 
duces more accurate ephemerides. 

• The accuracy of the ephemeris generated by the 
Chebyshev method increases with the degree of the 
Chebyshev polynomials very rapidly but bottoms 
out after a critical degree is reached. 

• The accuracy is generally better for smaller fit- 
ting spans. 

• The efficiency of the Chebyshev method is mainly 
related to the degree of the Chebyshev polynom- 
ials and the fitting span. 

• The Chebyshev method is slower than the Cowell 
method. Unless Chebyshev coefficients are re- 
quired, the Chebyshev method is not recommended 
for use in general applications. A study is cur- 
rently underway to further improve the efficiency 
of the Chebyshev method by using the Brouwer- 
Lyddane theory instead of the two-body theory for 
the starter. 

• A preliminary empirical formula was deduced to 
determine an optimal fitting span with a desir- 
able degree of Chebyshev polynomials in terms of 
high accuracy of the satellite ephemeris and low 
consumption of computer resources. 
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• It is also recommended that the conclusions of 
any study involving the use of Chebyshev ephem- 
erides obtained from the previous version of GTDS 
should be re-evalutated, especially in regard to 
accuracy. 
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APPENDIX A - MATHEMATICAL THEORY OF THE CHEBYSHEV 
ORBIT GENERATION METHOD 


A. 1 PROPERTIES OF CHEBYSHEV POLYNOMIALS 
The properties of the Chebyshev polynomials are most 
easily examined in the normalized interval [1, -1]. Any 
arbitrary finite interval [t^, t b l can be transformed 
to the normalized interval [1, -1] by the change of vari- 
able : ^ 



( A-l) 


where £ = the normalized time variable 

t = the start time of a polynomial fitting span, 
a (i.e., the start time of an integration step in 
GTDS terminology) 

t = the end time of a polynomial fitting span, 
b (i.e., the end time of an integration step and, 
therefore, t - t corresponds to the 
"step size"; a 

The Chebyshev polynomials are defined as a set of poly- 
nomials 


T^ (£) = cos i0 i = 0, 1, ... (A-2) 


1 The transformation could have been defined as 



so that t a would correspond to -1, and tb would corre- 
spond to +1. Since Reference 1 and the GTDS software have 
consistently used the definition as shown in Equa- 
tion (A-l) , this transformation is retained throughout 
this document and the new software. 
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generated from the sequence of cosine functions using the 
transformation 


6 = cos -1 ^ -1 £ £ £ 1 

Clearly for the zeroth degree 

T q (S) = cos (0) = 1 


(A-3) 


(A-4) 


and for the first degree 


T^d) = £ 


(A- 5) 


By repeated trigonometric manipulations, higher-degree 
Chebyshev polynomials can be computed yielding the recur- 
sion relation 

T.(£) = 2$T._ 1 U) - T i-2 (4) i = 2, 3, ... (A- 6) 

Table A-l contains the first ten Chebyshev polynomials. 

With simple algebraic manipulation, the algebraic func- 
tions, £ n , can be expressed in terms of a linear 
combination of the Chebyshev polynomials. This is shown 
in Table A-2. All the Chebyshev polynomials have a maxi- 
mum magnitude of 1 in the interval [1, -1] . The Chebyshev 
polynomials of degrees 0 to 3 are plotted in Figure A-l. 
The function of a parabola, | 2 , is also plotted in the 
figure as a linear combination of t q and T 2 . Except 
for Tq, all other Chebyshev polynomials cross the 
£-axis . The number of times that axis is crossed is equal 
to the degree of the Chebyshev polynomial and only those 
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Table A-l. The Chebyshev Polynomials 



T 2 = 2 £ 2 - 1 

t 3 = 4 i 3 — 3 a 

T 4 = 8 s 4 - 8 i 2 + 1 
T 5 = 16 ? 5 — 20 I 3 +5 ? 

Tg = 32 £ 6 -48 ? 4 + 18£ 2 - 1 
T 7 = 64? 7 -112{ 5 + 56? 3 -7? 

T d = 128 £ 8 -256? 6 + 160 f 4 - 32 t? + 1 

O 

T g = 256 f 9 - 576 $ 7 + 432 | 5 - 1 20 t 3 + 9 £ 
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Table A- 


. An Algebraic Function Expressed in Terms of a 
Linear Combination of the Chebyshev Ploynomial 


« - T i 

i 2 - (t 0 + t 2 )/2 

? 3 * (3T 1 + T 3 )/4 

J 4 = (3T 0 +4T 2 +T 4 )/8 

i 5 = (10T, +5T 3 +T 5 )/16 

£ 6 = (10T 0 + 15T 2 + 6T 4 + Tgl/32 

f 7 = (35T 1 + 21T 3 + 7T 5 + T 7 )/64 

£ 8 = (35T 0 + 56T 2 + 28T 4 + 8Tg + T g )/1 28 

t 9 = ( 1 26T, + 84T 3 + 36T 5 + 9T ? + T g ) /256 
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of odd degrees will cross the origin. The n locations at 
which the Chebyshev polynomial T n (£) crosses the l-axis 
are the n roots in the interval [1, -1] and are given by 


*k 


cos 


(2k - 1 ) 7T 
2n 


k = 1 , 


n 


(A-7) 


The significant properties of using the Chebyshev poly- 
nomials to fit an arbitrary function are that the error in 
the approximation is distributed evenly over the interval 
and the maximum error is reduced to the minimum or near- 
minimum value (References 3 and 4) . 

A. 2 INTERPOLATING POLYNOMIALS CONSISTING OF A LINEAR 
COMBINATION OF CHEBYSHEV POLYNOMIALS OF DIFFERENT 
DEGREES TO REPRESENT ACCELERATION 

Each component of the acceleration vector exerted on the 
spacecraft can be approximated by an interpolating poly- 
nomial consisting of a linear combination of Chebyshev 
polynomials: 


x(!) = p n m 


n 

E 

i=0 


C i T i 


( 4 ) 


(A-8) 


where x = the Cartesian component of the acceleration 
vector 

P n = the interpolating polynomial of degree n 

C- = Chebyshev coefficients for an acceleration 
component 

£ = the transformed time variable 

The accuracy of this approximation is better when the 
higher degrees of Chebyshev polynomials are included. 
However, the benefit of including higher degrees drops off 
quickly. This point is further illustrated in Section 4. 
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It is well known that the Chebyshev polynomials are one of 
the families which possess the property of orthogonality 
(References 3 and 4) . They are orthogonal in the interval 
[1, -1] with respect to the weighting function, 

w({) =1/ Jl - t 2 , i.e. , 



2 T i U) TjU) d| = 0, 
i-p |T.U)J 2 dl = A. t 0 


i + j 


(A-9) 


where A^ is a normalization factor which depends on i. 

Making use of the property of orthogonality, as demon- 
strated in Equation (A-9) , the Chebyshev coefficients can 
be evaluated 


C i 


1_ 

A. 




7^T 


i = 0, 1, 


T.(«) 


n 


xU) d£ 


( A-10) 


The above integral is difficult to evaluate because of the 
complexity of x(£) . However, it has been shown (Refer- 
ences 3, 4, and 5) that Equation (A-10) may be approxi- 
mated by 


C 


0 


1 

n + 1 


n+1 


E 

k=l 


x(4 k ) 


C. 

l 


2 

n + 1 


n+1 

E T i (5 k )x(€k ) 

k=l 1 k k 


i = 1, 2, 


(A-ll) 


n 
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where are the roots of the Chebyshev polynomial of 
degree n + 1, T n+1 (g). Therefore, with the accelera- 
tions evaluated at all the n + 1 roots, the variation of 
the acceleration in the interval [1, -1] , corresponding to 
the time interval [t a , t b ] , can be represented by the 
interpolating polynomial, P n (|), of degree n. 

A. 3 INTEGRATION OF CHEBYSHEV INTERPOLATING POLYNOMIAL TO 
. GENERATE EPHEMERIS 

With the acceleration components represented by Chebyshev 
interpolating polynomials as shown in Equation (A-8) , in- 
tegrating the equation once gives the velocity compo- 
nents, x: 


x({) = 


/ P " 


n 


(€) 


= £ 

i=0 


h 


(I) d$ 


(A-12) 


Through the use of Equations (A-4) , (A-5) , and (A-6) , the 
integration of the Chebyshev polynomials of different de- 
grees can be obtained: 


/v” 

/ T i <{) 


d4 = T 1 (4) + K 0 

d « - j[V*> + T 2 (i) ] + K 1 

dS = [i ttt T i+i ({) - r-h: T i-i <5) ] 

i — 2, 3, « • • , n 


(A-13) 


+ K. 
1 


t 


where K Q , K lf and are integration constants. 
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Substituting Equation (A-13) into Equation (A-12) and col- 
lecting terms of Chebyshev polynomials of the same degrees 
yields another interpolating polynomial of the following 
form for the velocity components: 


/ n+1 

P n l*> d* = Q n+i U) = Z) b i T i(0 (A-14) 


where Q , = the interpolating polynomial of degree n+1 


‘n+1 

b. 

1 


Chebyshev coefficients for a velocity compo- 
nent 


with 


b 0 K 0 + K i + ••• + 4 T o 


b l C 0 “ 2 C 2 


b . = ^-r 
i 2i 


C i-1 " C i+1 


(A-15) 


i = 2 , 3, . . . , n + 1 


C n+1 C n+2 0 


The integration constants in the expression for b Q may 
be evaluated from the initial velocity, i.e., x(t a ) = 
x( | = 1) : 

. n+1 

x( £ = 1) - E b t T.(| = 1) = b Q T Q U = 1) = b Q ( A-16) 

i=l 
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The interpolating polynomials for position components, 
x(£), can be obtained by the same procedure: 


- n+2 

*U> - _/o n+ l ( «> d « ■ R n + 2<«> - a l T i<«> 


(A-17) 


n+2 

= X( $ = 1) - £ a T U = 1) 

i=l 1 1 


a l b 0 " 2 b 2 


a . = 

l 


2i 


b i-l “ b i+l| 


i = 2, 3, . . . , n + 2 


b = b _ = 0 
n+2 n+3 


where R = the interpolating polynomial of degree n+2 
n L for the position component 

With velocity and position represented by Equations (A-14) 
and (A-17) in the interval [t t^] , the ephemerides 
of spacecraft at any other time within the interval can 
now be accurately and easily interpolated. 
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APPENDIX B ~ COMPARISON BETWEEN THE IMPROVED IMPLEMENTATION 
AND THE PREVIOUS IMPLEMENTATION OF THE CHEBYSHEV METHODS 


A series of GTDS computer runs was executed to compare the 
new and the previous software implementations in terms of 
their accuracy and efficiency. The accuracy was measured 
with respect to the ephemeris generated with the high- 
precision Cowell numerical integration method by using the 
GTDS Ephemeris Comparison Program. The efficiency is sim- 
ply a comparison of the CPU and I/O times consumed by the 
two different Chebyshev implementations. 

Three sets of test runs were made on GTDS with the GEOS-3 
satellite (arbitrarily chosen) over a two-day span using 
the Cowell method and the Chebyshev polynomial method of 
the new and previous implementations. The comparison 
results are presented in Table B-l. In these tests, the 
ephemeris generated by the Cowell integration method with 
a 24-second step size was used as a reference. The per- 
turbation (or force model) included in the Cowell method 
was identical to that used in the Chebyshev methods. 

Table B-l shows that the maximum difference in position 
vector of the ephemeris generated by the previous 
Chebyshev implementation with a 48th degree polynomial 
over a two-day arc is 97 meters with respect to the 
ephemeris generated by the Cowell method, while the new 
Chebyshev implementation with the same degree of poly- 
nomial has a maximum position difference of only 
0.25 meter. This represents an improvement of better than 
two orders of magnitude in the relative accuracy. 

The efficiency which is expressed as CPU time and I/O time 
consumed on the IBM S/360-75 computer is also examined. 
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Table B-l. Comparison Between New and Previous Chebyshev Implementations 
With Respect to the Cowell Method 
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COWELL ORBIT INTEGRATOR. 



The previous implementation used 8.012 minutes of CPU time 
and 0.154 minute of I/O time, while the new implementation 
used only 5.806 minutes of CPU time, a saving of 38 per- 
cent, and a comparable 0.142 minute of I/O time. 

The saving of computer resources can be viewed from 
another angle by lowering the degree of Chebyshev poly- 
nomials from 48 to 20 and 18. The results are also shown 
in Table B-l. Fitting Chebyshev polynomials with much 
lower degrees, the new implementation consumes four to 
five times less CPU resources yet maintains better ac- 
curacy than the previous implementation. 

Results in Table B-l indicate similar conclusions with 
more elaborate perturbation models, i.e., 8x8 geopotential 
field, 14th order resonance geopotential field, atmos- 
pheric drag, and solar radiation pressure as well as solar 
and lunar gravitational fields. 

The results presented in Table B-l are obtained with a 
fitting span of one satellite period for both the new and 
previous implementation. When the fitting span is in- 
creased to two satellite periods, the new implementation 
gives excellent results (I A R I = 0.12 meter). How- 
ever, after 80 loops in the iterative scheme, the previous 
implementation has simply failed to satisfy the 
10 -6 kilometer tolerance in fitting the first span and 
the computer run was subsequently terminated without gen- 
erating an ephemeris. 
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